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Ch : Abstract 

1—5 

^^ ■ We present a formahsm that improves the apphcabihty of perturbative QCD 
in the current region of semi-inclusive deep inelastic scattering. The formalism 

^^ , is based on all-order resummation of large logarithms arising in the pertur- 

(^ I bative treatment of hadron multiplicities and energy flows in this region. It 

OC ' is shown that the current region of semi-inclusive DIS is similar to the region 

\^ ' of small transverse momenta in vector boson production at hadron colliders. 

^D ■ We use this resummation formalism to describe transverse energy flows and 

0^ 



charged particle multiplicity measured at the electron-proton collider HERA. 
We find good agreement between our theoretical results and experimental data 



Q^' for the transverse energy flows. 



X 



I. INTRODUCTION 

It is well known that perturbative Quantum Chromodynaniics (pQCD) is a very powerful 
but not omnipotent theory of the strong interactions of elementary particles. It can be suc- 
cessfully applied to the calculation of various physical observables whenever the kinematics 
of the particle interaction process implies the existence of a large scale Q with dimension- 
ality of momentum. For such a kinematic regime, the cross-section of the hadronic process 
can be represented as a convolution of the perturbatively calculable hard part, describing 
the energetic short-range interactions of hadronic constituents (partons), and several process- 
independent nonperturbative functions, relevant to the complicated strong dynamics at large 
distances. 

The factorization of the hard and soft parts has proven to be a powerful method for the 
calculation of hadronic scattering cross-sections. Unfortunately, near the boundaries of the 
kinematic phase space the convergence of the perturbative solution can be spoiled by the 
presence of large logarithms logr, where r is some dimensionless function of the kinematic 
parameters of the system. For instance, r might be a small ratio of two momentum scales 
Pi and P2 of the system, r = P1/P2. 

To handle this situation, techniques for the all-order resummation of the logarithmically 
divergent terms have been developed |I|-^,|TB[. These techniques have been successfully used 
to improve the applicability of perturbative QCD in several processes (calculation of energy 
correlations in e~^e~ annihilation [^, transverse momentum distributions in vector boson 
Il5|j20|-p2| ,di-photon and Higgs production at hadron colliders). 



In this paper, we will consider another process, the production of hadrons in deep-inelastic 
lepton-hadron scattering (DIS). As will be discussed below, some features of this process 
are similar to the e~^e~ hadroproduction and vector boson production, so that in certain 
kinematic regions the description of this process requires all-order resummation of the large 
logs which would otherwise spoil the convergence of the perturbative calculation. 

Deep-inelastic lepton-hadron scattering (DIS) at large momentum transfer Q is one of 
the cornerstone processes to test pQCD. Traditionally, the experimental study of the fully 
inclusive DIS process, e + A -^ e + X where A is usually a nucleon, and X is any final 
state, is used to obtain information about the Parton Distribution Functions of the nucleon 
(PDFs). These functions describe the long-range dynamics of hadron interactions, and are 
required by many pQCD calculations. 

During the 1990 's, significant attention has also been paid to various aspects of semi- 
inclusive deep inelastic scattering (semi-inclusive DIS), for instance, the semi- inclusive pro- 
duction of hadrons and jets, e + A—^e + B + X and e + A —* e + jets + X. In particular, 
the HI and ZEUS collaborations at HERA, and the E665 experiment at Fermilab performed 
extensive experimental studies of the charged particle multiplicity [p HlO|] and hadronic trans- 
verse energy flows |TT[] at large momentum transfer Q. It was found that some aspects of the 



data, e.g., the Feynman x distributions, can be successfully explained in the framework of 
perturbative QCD analysis |T^. On the other hand, the applicability of pQCD for the de- 
scription of other features of the process is limited. For example, the perturbative calculation 
in lowest orders fails to describe the pseudorapidity or transverse momentum distributions 
of the final hadrons. Under certain kinematic conditions the whole perturbative expansion 



as a series in the QCD coupling may fail due to the large logarithms mentioned earlier. 

To be more specific, consider semi-inclusive DIS production of hadrons of a type B. At 
large energies, we can neglect the masses of the participating particles. In semi-inclusive DIS 
at given energies of the beams, any event can be characterized by two energy scales: the 
virtuality of the exchanged photon Q and the scale qt related to the transverse momentum 
of the final hadron B. The exact definition of qt will be given in the main part of the paper. 
One may try to use pQCD in any of three regions, where Q, qx, or both Q and qt are 
large. The renormalization and factorization scales should be chosen to be of the order of 
the large physical scale of the process. In the limit Q <^ Qt (photoproduction region) pQCD 
may fail due to the large terms (logQ/qT)"' as Q ^ 0, which should be resummed into the 
parton distribution function of the virtual photon |16|. The limit Q ^ qx is similar to the 



limit of a small transverse momentum in vector boson production where logarithms of the 
type (loggr/Q)" should be resummed in order to get a finite cross-section of the process 
1^. Finally, in the region qrp k. Q one may encounter another type of large logarithms 
corresponding to events with large rapidity gaps Ay. This type of large logarithms can be 



resummed with the help of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) formalism |]T3 . 

The purpose of this paper is to discuss the resummation of large logarithms in semi- 
inclusive DIS hadroproduction e + A^e + B + X'm the limit qx '^ Q- These are large logs 
arising from to the emission of soft and collinear partons, which we resum using the formalism 
of Collins, Soper, and Sterman (CSS) [^. Our calculations are based on the works of Meng, 



Olness, and Soper |jT^,|T5|], who analyzed the resummation technique for a particular energy 



distribution function of the semi-inclusive DIS process. This energy distribution function 
receives contributions from all possible final state hadrons, and does not depend on the 
specifics of fragmentation in the final state. 



In this paper we present a more general formalism than the one developed in [|I^,^, which 
will also account for the final state fragmentation of the partons. This formalism requires the 
knowledge of the fragmentation functions (FFs) describing nonperturbative fragmentation 
of final partons into observed hadrons. Correspondingly, this formalism can consider a wider 
class of physical observables including particle multiplicities. Our calculations will be done 
in the next-to- leading order of perturbative QCD. As an example of a practical application 
of our formalism, we compare our calculation with the HI data on the pseudorapidity distri- 



butions of the transverse energy flow |TT[ in the 7*p center-of-mass frame. We also present 
predictions for charged particle multiplicity. Another goal of this study is to find in which 
regions of kinematic parameters the CSS resummation formalism is sufficient to describe 
the existing data, and in which regions significant contributions from other hadroproduction 
mechanisms, such as the BFKL interactions [1T3, higher order corrections including multi- 



jet production with [jT6| or without ^7^ resolved photon contributions, or photoproduction 
showering |l^, cannot be ignored. 

The outline of the paper is as follows. In Section |I| we define the kinematic variables 
for the semi-inclusive DIS process and specify the coordinate frames which will be used 
throughout the following discussion. In Section ^ we derive the resummed cross-section 



formulas. In Section |^ we extend the results of Section |T| to obtain the resummed energy 



flows. In Section |V| we describe the matching between the resummed and perturbative cross- 
sections and energy flows. We also discuss kinematic corrections, which should be applied 
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FIG. 1. Semi-inclusive deep inelastic scattering 



to the resummed cross-section to account for fast contraction of the phase space of the 
perturbative cross-section at qt ~ Q- In Section |V| we describe the results of Monte Carlo 
calculations for the resummed cross-sections and energy flows. We present the comparison 
of our calculation with the existing energy flow data from HERA. We also suggest how one 
can reanalyze the existing HERA and Fermilab-E665 data on charged particle production in 
order to adapt it for unambiguous extraction of the non-perturbative Sudakov factors. 



II. KINEMATIC VARIABLES 



We follow notations which are similar to the ones used in IIT3,n!5i. In this Section we 



summarize them. 

We limit the discussion to the case of semi-inclusive DIS at the ep collider HERA. We 
consider the process 



e + A 



B + X, 
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where e is an electron or positron, A is a proton, i? is a hadron observed in the final state, 
and X represents any other particles in the final state in the sense of inclusive scattering 
(Fig. |I]). We denote the momenta of A and B by pj^ and p^, and the momenta of the 
lepton in the initial and final states by l^ and l''^. Also, g'^ is the momentum transfer to the 
hadron system, q^ = l^ — l'^ . For most of this paper, until the discussion of charged particle 
multiplicity, we will neglect the particle masses. 



We assume that the initial lepton and hadron interact only through a single photon 
exchange. Therefore, q^ also has the meaning of the 4-momentum of the exchanged virtual 
photon 7*. 

A. Lorentz scalars 

For further discussion, we define five Lorentz scalars relevant to the process (|ip. The first 
is the center-of-mass energy of the initial hadron and lepton ^JSeA where 

S,A={PA + lf = '^PA-l. (2) 

We will also use the conventional DIS variables x and Q^ which are defined from the mo- 
mentum transfer q^ by 

g2 = _q2 ^ 2£ . e, (3) 

X = ^ . (4) 

In principle, x and Q^ can be completely determined in an experimental event by measuring 
the momentum of the outgoing lepton. 

Next we define a scalar z related to the momentum of the final hadron state B by 

_ Pb -Pa _ 2xpB ■ Pa f^^ 

q-PA Q^ 

The variable z plays an important role in the description of fragmentation in the final state. 
In particular, in the quark-parton model (or in the leading order perturbative calculation) it 
is equal to the fraction of the fragmenting parton's momentum carried away by the observed 
hadron. 

The next relativistic invariant q^ is the square of the component of the virtual photon's 
4-momentum q^ that is transverse to the 4-momenta of the initial and final hadrons: 

4 = -^r^t/., (6) 

where 

I, a i, Q-Pb ^, q-PA .„x 

Qt =q^ - Pa Pb • iV 

Pa -Pb Pa- Pb 

The orthogonality of gf to both p^ and p'^, that is qt ■ Pa = qt ■ Pb = 0, follows immediately 
from its definition (^). 

The variable qx plays the same role in the semi-inclusive DIS resummation as the trans- 
verse momentum of a vector boson in resummation of vector boson production at hadron 
colliders. In particular, the theoretical cross-section calculated in a fixed order of pQCD is di- 
vergent in the limit gr —>■ 0, so that all-order resummation is needed to make the predictions 
of the theory finite in this limit. 




FIG. 2. Geometry of the particle momenta in the hadron frame 

In the analysis of kinematics, we will use two reference frames. The first is the center- 
of-mass (cm.) frame of the initial hadron and the virtual photon. The second is a special 
type of Breit frame which we will call, depending on whether the initial state is a hadron or 
a parton, the hadron or parton frame. As was shown in |]14|,|15[, by using the hadron frame 



one can organize the resummation formalism for semi-inclusive DIS in a way that is similar 
to the case of vector boson production. On the other hand, many experimental results are 
presented in the 7*p cm. frame. We will use a subscript h and cm to denote kinematic 
variables in the hadron or cm. frame. 



B. Hadron frame 



Following Meng et al. P| , p!5| the hadron frame is defined by two conditions: (a) the 4- 
momentum of the virtual photon is purely space-like, and (b) the momentum of the outgoing 
hadron B lies in the xz plane. The directions of particle momenta in this frame are shown 
in Fig. |. 

In this frame the proton A moves in the +z direction, while the momentum transfer q is 
in the —z direction, and g° is 0: 



g;: = (0,0,0,-Q), 
Pa,. = I (1,0, 0,1). 



(9) 



The momentum of the final-state hadron B is 

zQ f ql 2gT n ^T t\ n n^ 

P^B,. = ^(l + g„^,0,--l). (10) 

The incoming and outgoing lepton momenta in the hadron frame are defined in terms of 
variables ip and as follows |jl9|: 

£/^ = — (cosh?/', sinh-?/; cos 0, sinh-^ sin 0, —1), (11) 

£^ = — (cosh?/', sinh?/' cos ^jSinh?/' sin 0, +1). (12) 

Note that is the azimuthal angle of Ih or ^'h around the 2;-axis. -0 is a parameter of a boost 
which relates the hadron frame to a lepton Breit frame in which £^ = (Q/2, 0, 0, —Q/2). By 
(S) and (0) we find that 



cosh?/' = ^-l = --l, (13) 

Q^ y 

where the conventional DIS variable y is defined as 

y = -^- (14) 

The allowed range of the variable y in deep-inelastic scattering is < y < 1; therefore "0^0. 
The transverse part of the virtual photon momentum gf has a simple form in the hadron 
frame; it can be shown that 

<. = (-|,-gT,0,-|). (15) 

In other words, qt is the magnitude of the transverse component of qt^h- The transverse 
momentum pt of the final state hadron B in this frame is simply related to g^, by 

Pt = zqr- (16) 

Also, the pseudorapidity of B in the hadron frame is 

r]h = - log(tan -|^) = log -^. (17) 

The resummed cross-section will be derived using the hadron frame. To transform 
the result to other frames, it is useful to express the basis vectors of the hadron frame 
(T'^, X^, y^, Z^) in terms of the particle momenta [Q. For an arbitrary coordinate frame. 





Q ' 


X^ = 


qT\z 


Y^ = 


t^^'P'^ Z,TpX„ 


Z" = 


q. 



1+^ 

g2 



xp'X 



Q 



If these relations are evaluated in the hadron frame, the basis vectors T^, X^, Y^, Z^ are 
(1, 0, 0, 0), (0, 1, 0, 0), (0, 0, 1, 0), (0, 0, 0, 1), respectively. 

The relationships between the hadron-frame variables and the HERA lab-frame momenta 
are presented in Appendix A. 

C. Photon-hadron center-of-mass frame 



The center-of-mass frame of the proton A and virtual photon 7* is defined by the condition 
'PA,cm + Qcm = 0. The relationship between particle momenta in this frame is illustrated in 
Fig. ^. As in the hadron frame, the momenta qcm and pA.cm in the cm. frame are directed 
along the z axis. The coordinate transformation from the hadron frame into the 'j*p cm. 
frame consists of (a) a boost in the direction of the virtual photon and (b) inversion of the 
direction of the z axis, which is needed to make the definition of the cm. frame consistent 
with the one adopted in HERA experimental publications. In the cm. frame the momentum 
of 7* is 






w 



2W 



, u, u. 



2W 



(19) 



where W is the cm. energy of the 7*p collisions. 



w 



(pa + qf. 



(20) 




FIG. 3. Particle momenta in the 7*p cm. frame 



The momenta of the initial and final hadrons A and B are given by 

PA,cm-(^ 21^ '^'^' 2W )' ^ ' 

PB,cn. = (eb, PB,cm sin 0'^^ , 0, PB,c^ COS ^^ ] • (22) 

Since the hadron and cm. frames are related by a boost along the z-direction, the ex- 
pression for the transverse momentum of the final hadron B in the cm. frame is the same 
as the one in the hadron frame, 

Pt = zqt- (23) 

Also, similar to the case of the hadron frame, the relationship between qt and the pseudo- 
rapidity of B in the cm. frame is simple, 

q^ = We-''^"^. (24) 

The limit of small Qt, which is most relevant for our resummation calculation, corresponds 
to the region of large pseudorapidities in the hadronic cm. frame. Since in this case the 
final parton is produced in the direction of the momentum of the virtual photon, the region 
of large r]'^"^ is also called the current region. 

D. Parton kinematics 

The kinematic variables and momenta discussed so far are all laboratory variables. Next, 
we relate these to parton variables. 

Let a denote the parton in A that participates in the hard scattering, with momentum 
Pa — CaPA- Let b denote the parton of which i? is a fragment, with momentum p^^ = p%/^b- 
The momentum fractions ^a and ^h range from to 1. At the parton level, we introduce the 
Lorentz scalars £, z, qt analogous to the ones at the hadron level: 

X = = — , (25) 

2pa ■ q U 

. Pb-Pa Z 

z = = T^ (26) 

Q-Pa ^b 

Qt = -Qtqt,. (27) 

Here gj. is the component of q^ which is orthogonal to the parton 4-momenta p^ and p|^, 

qt ■ Pa = qt ■ Pb = 0. 
Therefore, 

-M u n q ' Pb a q ' Pa /no\ 

qt =q^ - p'i Pb • (28) 

Pa -Pb Pa- Pb 

In the case of massless initial and final hadrons the hadronic and partonic g^'s coincide, 

qT = qT- (29) 







T 





f) 



FIG. 4. Feynman diagrams for semi-inclusive DIS: a) LO; b-d)NLO virtual diagrams; e-f) 
NLO real emission diagrams 
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III. THE RESUMMED NLO CROSS-SECTION 

The knowledge of five Lorentz scalars SeA, Q, Qt, x, z and the lepton azimuthal angle 
(j) in the hadron frame is sufficient to specify unambiguously the kinematics of the semi- 
inclusive scattering event e + A~^e + B + X. In the following, we will discuss the hadron 
cross-section das a, which is related to the parton cross-section data by 

dxdzdQMqld<P ^Jz i, '"'^^"'^'''J. ^a ^^^^^'^^'dxdzdQMq^df ^ ' 

Here Fa/A{C,a, fJ^r) denotes the distribution function (PDF) of the parton of a type a in the 
hadron A, and Ds/bi^b, f^'o) is the fragmentation function (FF) for parton type b and final 
hadron B. The parameters fip and fio are the factorization scales for the PDFs and FFs. 
In the following discussion and calculations, we assume that these factorization scales and 
the renormalization scale fiR are the same. 

The analysis of semi-inclusive DIS can be conveniently organized by separating the de- 
pendence of the parton and hadron cross-sections on the leptonic angle and the boost 
parameter ip from the other kinematic variables x, z, Q and qt [19|. Following |n], we 



express the hadron (or parton) cross-section as a sum over products of functions of these 
lepton angles in the hadron frame y4a(?/^,0), and structure functions Vg*^ (x, 2;, Q^, g|.) (or 
yba\^, z, Q^ ^t), respectively): 



dasA 



dxdzdQ'^dqj^dcj) 



Y^vl,i{x,z,Q\ql)A^{^,4>), (31) 



a=l 



d(Tba _ sr^ f^(a) 



dxdzdQ'^dq^dip 



^n^(x,^,Q^g|)A„(^,</)). (32) 



a=l 



At the energy of HERA, hadropro duct ion via parity-violating Z-boson exchanges can be 
neglected, and only four out of the nine angular functions listed in [|I^ contribute to the 
cross-sections (^l]-^2|). They are 



v4i = 


1 + cosh ■?/', 


A2 = 


-2, 


A, = 


— cos (f) sinh 2il) 


A4 = 


cos20sinh ip. 



(33) 

We will assume that the angle is not monitored in the experiment, so that it will 
be integrated out in the following discussion. Correspondingly, our numerical result for 
da / {dxdzdQ'^dq^) will not depend on terms in (|3l| - |32D proportional to the angular functions 
^3 and A4, which integrate to zero. 
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Out of the four structure functions, only V^^ receives contributions from both leading 
and next-to-leading order diagrams. Also, only the V^^ structure function diverges in the 
limit qt ^ 0- 

The leading order (LO) parton process is e + a^e + b where the initial parton a is a 
quark from the proton and the final quark b (which is the same as a) fragments into the 
hadron B. The Feynman diagram for this process is shown in Fig. ^. There is no LO 
contribution to semi-inclusive DIS from gluons. 

The next-to-leading order (NLO) corrections are shown in Figs. ^-f. At this order, we 

need to account for the virtual corrections to the LO subprocess Q Y ^ ^ (Figs, ^-d), as 

well as for the diagrams describing the subprocesses Q 'y* —>■ q g and g'-f* —>■ qq, with the 
subsequent fragmentation of the final-state quark, antiquark or gluon (Figs. ^-f). 

Conservation of total 4-momentum in the real emission subprocesses (Figs. §e-f ) allows 
us to write the momentum of the unobserved final state parton (e.g. the gluon in Fig. ^) as 

P'g=q'' + P'a-Pb- (34) 

When there is no gluon radiation (p^^ = 0) the momentum of bis p^^ = Pa+q'^, and, according 
to (|28|) , q^ = —qt ■ qt = 0. Thus, a non-zero qx in the event is an effect of gluon radiation. 
In the region qx —>■ 0, either softness or coUinearity of the unobserved partons will create 
infrared singularities, which make the perturbative result unreliable. The sum of the real 
and virtual diagrams is made finite by order-by-order cancellation of the soft singularities 
arising from the real and virtual pieces, and by absorption of the collinear singularities into 
the parton distribution and fragmentation functions. Nonetheless, this cancellation does 
not guarantee rapid convergence of the perturbative calculation, which will typically contain 
large logarithms log qr/Q countering the smallness of the strong coupling. 

The slow convergence of the perturbative series at gr — ^ can be corrected by resum- 
mation of the most singular logarithmic terms. It is done in the following way. First, we 
extract the terms in the squared amplitudes of the real emission diagrams Figs, ^e-f that are 
most singular in the limit g-j- ~^ 0; we refer to these terms as the asymptotic piece. These 
terms are proportional to l/q^ and, as it was mentioned above, they appear only in the V^^ 
structure function. Thus, the structure function V^J is represented as 



V,':>ix,z,Q',q'r) = (V,':\x,z,Q,q'j.) + Y^^\x ^z . Q\ q^) . (35) 

\ / asyva 

where {V^atasym is C(l/g|'), and 1^^ is finite in the limit qx -^ 0. The asymptotic piece of 
the NLO hadron cross-section (|30|) is 



\dxdzdQ'^dq^d(f) J SeA tt 2q^ 2% 



asyva 



E-l 



.2 

i 



DB/j{Z,^l)[{Pgg O fj/A){X,fi) + {Pgg O f g / a) {x , fl)} 



+ {{Db/j O Pgg){z,n) + {Db/s O Pgg){z,fi)}fj/A{x,fi) 



+ 2DB/A^,fi)f,/A{x,fi){CFlog%-^Cp} + 0{i^^) ,ql 
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(36) 



Here ecj is the electric charge of the participating quark or antiquark of flavor j; the param- 
eter o"o is 

Q2 .e^ 



the factor F;, that comes from the leptonic side, is defined by 

F, = ^^; (38) 

the color factor Cp = {N^ — l)/(2A^c) = 4/3. The convolution in (|36|) is defined as 

(/o^)(a;,/x) = £/(x/e,/ik(e,/i)y • (39) 

The functions Pij{x) entering the convolution integrals in (|36D are the familiar splitting 
kernels: 

P,,{x)=Cp(^^^ , (40) 

P,,ix) = ^(1 - 2x + 2x''^ , (41) 

P,,ix) = C^l±ili:^. (42) 

X 

The finite piece YgJ of the hadron cross-section and the other structure functions V^„ for i = 
2,3,4 can be derived in a straightforward way from the expression for the cross-section of 
the real emission subprocesses, which is presented in Appendix B. 

Next, we use the perturbative asymptotic piece (|36| ) to derive the 0{as) expression 
for the resummed cross-section {da / {dxdzdQ'^dq^d(f)))resum- In the CoUins-Soper-Sterman 
resummation formalism f^, the cross-section is written as a Fourier integral over a variable 
b conjugate to qt, 

dasA aoFiA,iij,(P) r d^-% 



/do 
7^:;T^e''^-'WBA{b, X, z, Q) + Yba- (43) 



dxdzdQ"^ doj-dip S^a 2 J (2-7? 

The term containing the integral of WBA{b), which we name the CSS piece, absorbs the 
asymptotic contributions from all orders. The second term, which is the finite piece, has the 
form 



Yba = 4'i + E Vs'iM^, 0)- (44) 



The Fourier transform is performed in the space of n = 4 — e dimensions, in which the 
asymptotic piece of the real emission subprocesses generates terms which are proportional 
to 1/e and divergent as e ^ 0. Upon summation of the real and virtual diagrams, some of 
these terms, specifically those corresponding to soft singularities, cancel. The remaining 1/e 

13 



poles correspond to collinear singularities. They are later absorbed into the redefined NLO 
PDF or FF, rendering a final expression that is nonsingular as e — *> 0. 

According to 0, the form of the resummed structure function WBA{b,x,z,Q) at small 
values of b is 

j 

In the limit of small b and large Q, the Sudakov function SsA^b, Q) does not depend on the 
types of the external hadrons and looks like 

rCiQ^ r/772 / C^O^ \ 

SsAib, Q)= I ^ Aiasi-p), Ci) In -^ + i?(a,(7l), C,, C,) , (46) 

with 

A(«.(7I), = f A,(Ci) (^^ , (47) 

B{aM, Cu C2) = E B,{C„ C2) (^) . (48) 



The integration in (^q) is performed between two scales of order 1/6 and Q. The constants 
Ci and C2 determining the exact integration range are a priori unknown, and their variation 
allows us to test the scale invariance of the resummed cross-section. There are convenient 
choices of Ci, C2 for which some logarithms in Wba{b, x, z, Q) cancel. 

The functions C^'"{x,b, n) and C°'^'^{z,b, fi) contain contributions from partons radiated 
collinearily to the initial and final hadrons. These functions can also be expanded in series 
of Os/tt, as 

Q;(£,&,/i) = f:Clf\S:,f,b)(^^]\ (49) 

k=o \ ^ J 

k=0 \ ^ / 

The renormalization scale in the C-functions is 

/i = 2e-yb = bo/b, 

where 7 = 0.577215... is the Euler constant. 

The explicit expressions for Ak{Ci), BkiCi, C2) and the C-functions can be obtained by 
comparing the expansion of WsAib^x^z.Q) as a series in Os/vr with the 6-space expression 
for the perturbative cross-section. Using our NLO results, we find 

Ai = Cp, (51) 

To the same order, our expressions for the C-functions are 
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LO: 



C;f\x,iih) = 5,k5{l-x), (53) 

C;f'\z,^^h) = 6,,5{l-z), (54) 



NLO: 



c;f^(£,H = ^(i-£)-i^.,(x)iog(^) 



C^5(l-x)(^ + log^(^-^)), (56) 



,16 "V boC2 
C;f)(x,/i6) = lx(l-x)-P,G(x)log(^), (57) 

-C..(l-.)(|.log^(^)), (58) 

C;rH^,/^&) = §^^-^G.(^)log(^). (59) 

In these formulas, the indices j and k correspond to quarks and antiquarks, and g to gluons. 
Due to the crossing relations between parton-level semi- inclusive DIS, vector boson produc- 
tion, and e^e^ hadropro duct ion, the C™-functions are essentially the same in semi-inclusive 
DIS and the Drell-Yan process; and the C°"*-functions are essentially the same in semi- 
inclusive DIS and e~^e~ hadropro duct ion. At NLO the only difference stems from the fact 
that the momentum transfer g^ is spacelike in DIS and timelike in the other two processes. 
The virtual diagrams Figs. ^D-d differ by tt^ for spacelike and timelike g^. Correspondingly, 
Cjl and C°^ for semi-inclusive DIS ( |5B| ) do not contain the term {tt'^/3)5{1 — x) (or 

(7r^/3)5(l — z) ) which is present in the C*^ -function for the Drell-Yan process (or in the 
C°^ -function for e^e" hadroproduction, respectively )Q On the other hand, the NLO ex- 
pression for the Sudakov factor (^6[) is the same for semi- inclusive DIS, the Drell-Yan process, 
and e^e~ hadroproduction, which also results from the crossing symmetry. 

Up to now, we have been discussing the behavior of the resummed cross-section at short 
distances. The representation (EH) should be modified at large values of the the variable h 



^With two minor exceptions, our expressions for the functions C are equivalent to the ones pub- 
lished previously by Meng et al. |jl5[. The 7r^/3 terms are incorrectly included in equations (43) 



and (45) of ||T^ for the semi-inclusive DIS functions C*^ and C°^ . Also, our Eq. ( [5^ ) contains 



23Ci7'/16 = 23/12 instead of 29/12 in Eq. (45) of [15|, which is apparently due to a typo in |K 
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to account for nonperturbative long-distance dynamics. The modified ansatz for Wba valid 
at all values of b is 



WBA{b,x,z,Q)=Y.e%DB/boC",f){z,K,fi){q-,oFa/A){x,K,fi) e- 

3 

Here the variable 

h 
6* 



jl I ( b V' 

y \ Omax / 



(60) 



(61) 



serves to turn off the perturbative dynamics for b > b^ax, with b^ax ~ 1 GeV~^. Fur- 
thermore, the Sudakov factor is modified, being written as the sum of the perturbatively 
calculable part 5''^(6^,,Q) given by (PB|), and a nonperturbative part which is only partially 
constrained by the theory: 



SsAib, X, z,Q) = S (6*, X, z, Q) + 5*^^ (6, x, z, Q). 



(62) 



From the renormalization properties of the theory, it can be concluded that the Q dependence 
of the nonperturbative Sudakov term should be separated from the dependence on the other 
kinematic variables, i.e. 



S^^{b,x,z,Q) 



9BAib, X, z) + g'i\{b) log — 



(63) 



with Qo ~ 1 GeV. The theory does not predict the functional forms of g^j^{b,x, z) and 
5'ba(^)' so these must be determined by fitting experimental data. Also, 5*^^ can depend on 
the types of the hadrons A and B. On the other hand, due to the crossing symmetry between 
semi-inclusive DIS, the Drell-Yan process and e+e" hadropro duct ion, one may expect that 
the functions g^'^\b) in these processes are not independent ||15|, but satisfy 



9B]i{b)l^,, = l{9'Akb)ly + Ab) 



(64) 



If the relationship (|6^ is true, then the function gBA^b) in semi- inclusive DIS is completely 



(2) 



(2) 



.(2) 



known once the parameterizations for the functions g^jiib) in the Drell-Yan and gBB^b) in 
e^e~ hadropro duct ion processes are available. In practice, the Drell-Yan nonperturbative 



Sudakov factor is known only when the incoming particle is a nucleon P0|-Pi| , while the non- 
perturbative Sudakov factor for e^e~ hadropro duct ion is available only for energy correla- 
tions IQ . An additional complication comes from the fact that the known parameterizations 
of the nonperturbative Sudakov factors for the Drell-Yan |2^-^ and e^e~ hadroproduction 



24 1 processes correspond to slightly different scale choices. 



Ci — ^07 C2 



and 



Ci 



Co 



-3/4 



(65) 



(66) 
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respectively. Therefore, the known functions g^'^' (b) and g^'^' _ (b) are not 100% com- 
patible, and in principle should not be combined to obtain g^'^\b) for semi-inclusive DIS. 

Despite this minor incompatibility, we will use (Q) to construct g^^\b) for our numerical 
calculation of energy flows and, with less justification, particle multiplicities. We have found 
that the numerical results for the energy flows are only slightly dependent on the choice be- 
tween the two sets (^^Hj) of the constants Ci, C2 (see Section ^I|). Also, detailed information 



about the functional form of the Q-dependent part of the nonperturbative Sudakov factor can 
be obtained only by studying the dependence of the resummed cross-sections on Q. Since the 
HERA data discussed in this paper covers only a small range of Q (3.62 < Q < 5.71 GeV), 
it is hard to distinguish between the uncertainties in the Q-dependent and constant parts 
of the nonperturbative Sudakov factor. Of course, more definite conclusions about g^'^^ will 
be possible once more detailed semi-inclusive DIS data from HERA and Fermilab-E665, 
covering a wider range of Q, are available. 

IV. HADRONIC MULTIPLICITIES AND ENERGY FLOWS 

Knowing the hadron cross-section, it is possible to calculate the multiplicity of the process, 
which is defined as the ratio of this cross-section, and the total inclusive DIS cross-section 
for the given leptonic cuts: 



1 da 

datot/dxdQ"^ dxdzdQ'^dq^ 



Multiplicity = ^. /j^j^, j^j..j^,.j 2 - (67) 



Both the cross-section and the multiplicity depend on the properties of the final-state frag- 
mentation. The analysis can be simplified by considering energy flows which do not have 
such dependence. A traditional variable used in the experimental literature is a transverse 
energy flow (-Et) in one of the coordinate frames, defined as 

{Et)^s = E/ d<^BET^ — '-. (68) 

This definition involves an integration over the available phase space <I>b and a summation 
over all possible species of the final hadrons B. Since the integration over $b includes 
integration over the longitudinal component of the momentum of -B, the dependence of {Ex) 
on the fragmentation functions drops out due to the normalization condition 

E / ^ DB/b{z)dz = 1. (69) 

B •' 

Instead of {Et)-, we analyze the flow of the variable z. This flow is defined as 

d^z V- f^ Me + A^e + B + X) 

-22 ^ JT^TTT^T^O::! dz. (70) 



dx dQ"^ dq^ ^ Jzmtn dx dz dQ'^ dq^ 
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We prefer to use S^ rather than {Et) because {Et) is not Lorentz invariant, which com- 
phcates its usage in the theoretical analysis^. Also, the analysis in terms of qt and 2;-flow 
makes the analogy between resummation in the current region of semi-inclusive DIS and in 
the small transverse momentum region of the Drell-Yan process more obvious. Since qt is 
simply related to the pseudorapidity in the 7*p cm. frame via Eq. (0), and the transverse 
energy of a nearly massless particle in this frame is given by 



^7 



(71) 



the experimental information on dTjz/{dx dQ"^ dq^) can be derived from the 7*p cm. frame 
pseudorapidity (rjcm) distributions of {Et) in bins of x and Q^. If mass effects are neglected, 
we have 



d{ET) 



Qt 



2 ULj; 



dricm dqr 

The asymptotic contribution to the 2;-flow distribution is 

f dE, \ ^ aoFias e^ Ai{^,(j)) 

asym 



dxdQ'^dqj^ 



SeA TT 2g|, 2-71 



Ee? 



{(Pqq ° fj/A){x, fi) + {Pgg O fg/A){x, /i)} 



+ 2/,/^(a;,/x){C^log^-|c^} + 0( 



a. 



vr 



Qt 



The resummed z-flow distribution is 

t/S, aoFiAi{ij,(P) r d^-% 



dxdQ'^dqj.i 



S. 



(2vr; 



n-2 



e"^'^-'W,{b,x,Q)+Y,. 



with 



W^ib, x,Q) = J2 e]CT\K,^i){C% o F,/^)(x, 6,, /i) e 



-S46,x,Q) 



(72) 



(73) 



(74) 



(75) 



In (I75|), the NLO function C°''\h,^i) is 



CT\h,^^) 



a. 



n 



^ / 65 4, fe/i 
^^(-48 + 3^°< 



log' 



C260 



)' 


asCp 


H 

vr 2 



1 8 , 6/i 

o - Q log IT 
3 3 Oq 



Similar to (0), the 2;-flow Sudakov factor 5*^ is a sum of perturbative and nonperturbative 
parts. 



5,(6, X, Q) = 5^(6,, X, Q) + 5f ^(6, x, Q) 



(76) 



^The z-flow E, is related to the energy distribution function S calculated in |15| as S, 
{2xEa/Q'^)'^- Here Ea is the energy of the initial hadron in the HERA lab frame. 



The NLO perturbative Sudakov factor 5*^ is given by the universal x-independent expression 
(^). By the same argument as in the case of semi- inclusive DIS multiplicities, we assume 
that the non-perturbative part of Sz can be parameterized as 

Sr{b,x,Q) = 9^'Hb,x) + l{9^'\b)l^ + 9^'\b)l^J log^. (77) 

In the numerical calculation, we use the functions g^'^\h) from [2C] and g^'^\h) _ from 



2^ , despite the fact that g^'^\b) was fitted to Drell-Yan data using a different C2 value 



DY 



than g^'^\b) 

We also parameterize the functional form of g^^\b, x) in terms of b and x as 

g^'\b,x) = ih + ^)b^ (7^ 

\/x 



where the constants hi, /12 must be determined by fitting the experimental data. 

In principle, the 2;- flow Sudakov factor Sz{b,x,Q) is related to the Sudakov factors 
SsAip, X, z, Q) of the contributing hadroproduction processes e + A^i-e + B + X through 
the relationship 



e-^^^"'^^ = ^^^ -Y. I ^dze-'''''^''^''''^\DB/boC",f){z,b,,fi). (79) 



-S,ib,x) ^ _ 



In practice, the efficient usage of this relationship to constrain the Sudakov factors is only 
possible if the fragmentation functions and the hadronic contents of the final state are accu- 
rately known. We do not use the relationship ( |79D in our calculations. 

V. RELATIONSHIP BETWEEN THE PERTURBATIVE AND RESUMMED 
CROSS-SECTIONS. UNCERTAINTIES OF THE CALCULATION 

In the numerical calculations, some care is needed to treat the uncertainties in the defini- 
tions of the asymptotic and resummed cross-sections (|36D and (|43|) , although formally these 
uncertainties are of order 0((a;s/7r)^, g^^). 



A. Matching 

The generic structure of the resummed cross-section (^), calculated up to the order 
0{{ajnr), is 



a. 



■Zim = ^css + Y^''\ (80) 



In (|80|), the CSS piece receives all-order contributions from large logarithmic terms 

„n 2n-l 



a 



-fY.^n.nlog'^iq'T/Q'), n=l,...,oo. (81) 

^T m=0 
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FIG. 5. 
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a) b) 

Choice of the matching between the resummed and perturbative cross-section 



The y-piece is the difference of the fixed-order perturbative and asymptotic cross-sections, 



y(A^) 



'-'•pert ^ asym- 



12) 



In the small-gr region, we expect cancellation up to terms of order 0{a^^^ /tt^^^) be- 
tween the perturbative and asymptotic pieces in (^2]), so that the CSS piece dominates the 
resummed cross-section (^Of). On the other hand, the expression for the asymptotic piece 
coincides with the expansion of the CSS piece up to the order 0{a^ /ti^). Therefore, at 
larger g-p, we expect better cancellation between the CSS and asymptotic pieces, and at 
large qt the resummed cross-section (^3|) should be equal to the perturbative cross-section 
up to corrections of order 0{a^^^ /t[^^^). 

In principle, due to the cancellation between the perturbative and asymptotic pieces at 
small qti and between the resummed and asymptotic piece at large qt-, the resummed formula 
Cresnm is at Icast as good an approximation of the physical cross-section as the perturbative 
cross-section apert of the same order. However, in the NLO calculation at large qt it is safer 
to use the fixed order cross- section ( ^OD instead the resummed expression (|^). At the NLO 
order of a^, the difference between the CSS and the asymptotic pieces at large g-r may still 
be non- negligible. Therefore, the resummed cross-section (Tresum niay differ significantly from 
the NLO perturbative cross-section (Tpert- This difference does not mean that the resummed 
cross-section agrees with the data better than the fixed-order one. At qt ~ Q-, the NLO 
cross-section is no longer dominated by the logarithms that are resummed in (|43|) . In other 
words, the resummed cross-section (^) does not include some terms in the NLO cross-section 
that become important at qt ~ Q- For this reason, at large qt the resummed cross-section 
may show unphysical behavior; for example, it can become negative or oscillate. 

As the order of the perturbative calculation increases, we expect the agreement between 
the resummed and the fixed-order perturbative cross-sections to improve. Indeed, such 
improvement was shown in the case of vector boson production [Q , where one observes a 
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smoother transition from the resummed to the fixed-order perturbative cross-section if the 
calculation is done at the next-to-next-to-leading order. Also, at the NNLO the switching 
occurs at larger values of the transverse momentum of the vector boson than in the case of 
the NLO. 

The switching from the resummed to the fixed-order perturbative cross-section should 
occur at qt ~ Q- Nonetheless, there is no unique prescription for the exact point at which 
it should happen. In our program we switch from the resummed cross-section to the pertur- 
bative result at the first minimum of the difference 

I da \ I da \ 



\ "^ / pert 

lying at qt above the main maximum of the resummed cross-section (Fig. ^. This prescrip- 
tion is satisfactory for two different possible situations, in which the resummed cross-section 
crosses (Fig. ^ja) or does not cross (Fig. ^d) the perturbative cross-section. 

B. Kinematic corrections at qt ^ Q 

In this subsection we will discuss the differences between the kinematics implemented 
in the definitions of the asymptotic and resummed cross-sections (|36| ) and (^), and the 
kinematics of the perturbative piece at non-zero values of qt- 

Consider the perturbative hadronic cross-section (^). The integrand of ( pUD contains 
the delta-function 






xz6 



i^a-x){^b- Z) -XZ-^ 



^3) 



which comes from the parton-level cross-section ( p.O(j| ). Depending on the values of 
X, z, Q^, g|^, the contour of the integration over ^a and ^b determined by (|83D can have one of 
three shapes shown in Fig. Pa,b,c. For qr "^ Q the integration proceeds along the contour 
in Fig. ^a, and the integral in ( pO]) can be written in either of two alternative forms 

das A [^ d^a 



dxdzdQ'^dq^dcp J{ia)rmn ia-x 



/ _" MBA{^a,^b;X,Z,Q^,qT, 



1 df, 

-^Mba{U 6; X, z, Q\ ql, 0), (84) 

where 

MsA^ia, ib] X, z, Q^ q^, 0) = 

j^f^,-xz Y: epB/b{^b)F./A{Q E f!;:\£, z, Q\ g|)A.(^, 0). (85) 

The lower bounds of the integrals are 
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w'^ 



{.ih)min = \- Z, (87) 

1 — X 



with 



w = ^yi^. (88) 

Alternatively, the cross-section can be written in a form symmetric with respect to x and 

z, 

dxdzdQ'^dq^d(f) Jx+w ia — x 

+ r -^M(ea,e6;x,5,Q',g|,0), (89) 

J z+w C,b — Z 

where the integrals are calculated along the branches RP and RQ in Fig. ^, respectively. 

As qt -^ 0, 



\<,a)min ^ -^i \Si 



bjmin 



and the contour PRQ approaches the contour of integration of the asymptotic cross-section 
(p6|) shown in Fig. ^. The horizontal (or vertical) branch contributes to the convolutions 
with splitting functions in (p6D arising from the initial (or final) state collinear singularities, 
while the soft singularities of (|36|) are located at the point ^a = x, ^b = z. 

On the other hand, as qt increases up to values around Q, the difference between the 
contours of integration of the perturbative and asymptotic cross-sections may become signif- 



icant. First, as can be seen from (|89|), in the perturbative piece ^a and ^b are always higher 
than X + w or z + w , while in the asymptotic piece they vary between x or z and unity. At 
small X (or small z) the difference between the phase spaces of the perturbative and asymp- 
totic pieces may become important due to the steep rise of the PDFs and FFs in this region. 
Indeed, for illustration consider a semi-inclusive DIS experiment at small x. Let qr/Q = 0.5, 
z = 0.1, and x = 10~^; then x + w = 1.6 ■ 10~^ 3> a; = 10~^. In combination with the fast 
rise of the PDFs at small x, this will enhance the difference between the perturbative and 
asymptotic cross-sections. 

Second, for x or z near unity, it could happen that x + w>1oyz + w>1, which would 
lead to the disappearance of one or two branches of the integration of the perturbative piece 
(Fig. P3,c). In this situation the phase space for nearly collinear radiation along the direction 
of the initial or final parton is suppressed. Again, this may degrade the consistency between 
the perturbative and asymptotic piece, since the latter includes contributions from both 
branches of the collinear radiation. Fortunately, the x — z asymmetry of the phase space in 
semi-inclusive DIS is not important in the analysis of the existing data from HERA, since it 
covers the small- a; region and is less sensitive to the contributions from the large z region, 
where the rate of the hadroproduction is small. However, in the numerical analysis we found 
it necessary to correct for the contraction of the perturbative phase space described in the 
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previous paragraph. We incorporate this correction by substituting for x and z in (^, ^3]) 
the rescaled variables 

. = -^^^- (90) 

These substitutions simulate the phase space contraction of the perturbative piece. At small 
qr, the rescaling reproduces the exact asymptotic and resummed pieces (^) and (^31), but 
at larger q^ it excludes the unphysical integration regions of .^a ~ a; and C,b ^ z. 
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FIG. 6. The contours of the integration over ^^j ?fe for a,b,c) the perturbative cross-section 
; d) the asymptotic and resummed cross-sections (36) and 
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VI. NUMERICAL RESULTS 

In this section, we present the results of Monte-Carlo simulations for the z-flow and the 
differential multiplicity of charged particle production, 



1 



t/S, 



and 



^^chgd 



datot/idxdQ"^) dxdQ'^dqr datot/{dxdQ'^) dxdzdQ'^dqx 

Our calculations use the parameters of the HERA electron-proton collider. The energies of 
the proton and electron beams are taken to be equal to 820 and 27.5 GeV, respectively. 



A. Energy flows 



As a first application of the resummation formalism, we consider the cm. pseudorapidity 
distributions of the transverse energy flows in the current region, data for which has been 
published in |11]]. We discuss the data in seven bins of x and Q, four of them covering the 
region 10 < Q^ < 20 GeV^ 3.7 ■ 10"^ < x < 2.3 ■ 10-^ and the other three the region 
20 < Q2 < 50 GeV^ 9.3 ■ lO"^ < x < 4.9 ■ lO^^. 



CTEQ4M 

X = 0.0049 
Q = 5.7 GeV 




FIG. 7. Comparison of the NLO perturbative (111) and resummed ( [7^ ) expressions for the 
z-flow distribution with the existing experimental data from HERA |T^ . The presented data is for 
(x) = 0.0049, (Q2) = 32.6 GeV^. 
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FIG. 8. Comparison of the resummed (|74| ) z-flow distribution with the HERA data from |11| 
in 7 bins of (x) and (Q). The resummed z-flows were calculated using the parameterization (pl|) 
for the non-perturbative part of the Sudakov factor. 



In our calculation, we use the CTEQ4M parton distribution functions |^. The factor- 
ization and renormalization scales of the perturbative and asymptotic pieces are all set equal 
to ^ = Q. As was mentioned in Section |I^, the data on transverse energy flow can be easily 
transformed into the qt distributions of the z-fiow. In Fig. |^, we present the comparison of 
the existing data in one of the bins of x and Q"^ {{x) = 0.0049, (Q^) = 32.6 GeV^) with 
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the NLO perturbative and resumined z-flows given in Eqs. (|111|) and (|74D, respectively. In 



Fig. ^ we present the comparison of the resummed 2;-flow with the data in the other bins of 



11 



Figure ^ demonstrates two important aspects of the NLO qr distribution, namely, that 
the NLO exceeds the data at small qt and is below the data at qx > Q. In fact, we find 
that the deficit of the NLO prediction of perturbative theory in comparison with the data 
at medium and large qx {qr > 5 GeV) is present in the entire region of x and Q^ that we 
have studied. 

As we discussed in Section |V], one can trust the resummed calculation only for reasonably 
small values of qr/Q- For large values of qx, the fixed-order perturbative result is more reli- 
able. This means that the NLO resummation formalism will not give an accurate description 
of the data for qx ^ Q, due to the small magnitude of the NLO perturbative 2;-flow in this 
region. 

The excess of the data over the NLO calculation can be interpreted as a signature of 
other intensive hadroproduction mechanisms at cm. pseudorapidities rj'^"^ < 2. A discussion 
of the cross-sections in this pseudorapidity region is beyond the scope of our paper. However, 
we would like to point out that there exist several possible explanations of the data in this 
region, for instance, the enhancement of the cross-section due to BFKL showering |TB[ or 
resolved photon contributions |16;18]. From the point of view of our study, it is clear that 
good agreement between the data and the combination of the perturbative calculation and 
the CSS resummation, in a wider range of t]*^™", will be achieved when next-to-next-to- leading 
order contributions, like the ones contributing to (2+1) jet production |l^, are taken into 
account. 

On the other hand. Figs. ^ and ^ illustrate that the resummation formalism accurately 
describes the data in the region qx < 10 GeV. This calculation of the resummed 2;-fiows ( |7^ 
was done using the following parameterization for the non-perturbative Sudakov factor (|77D 

S^''{b,x,Q)=g^'\b,x)+g^'\b,Q), (91) 

where the Q-dependent part g^^\b, Q) is completely defined by the symmetry between semi- 
inclusive DIS, Drell-Yan and e~^e~ hadroproduction processes (Section |T|), 



g^'\b, Q) = \b^ 0.48 log(-^) + 5.32C^ log(^) log(-^) . (92) 



The parameterizations of the Q-dependent parts of the non-perturbative Sudakov factors in 
the Drell-Yan and e^e~ hadroproduction processes are taken from [^ and |2^, respectively. 
In (p^), the constants are Ci = 2e~"', C2 = e~^^^,Qo = 1 GeV. The variable 6^, is given by 
(0), with bmax = 1 GeV-^ 

The x-dependent function g^^\x,Q) was determined by fitting the HERA data of [jll]. 



We found that good agreement with the data is obtained when g^^\b, x) is parameterized by 
a linear function of x~°'^, 

^(1) (b, x) = 6" f-2.58 + ^) , for X < 10'^^ (93) 
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FIG. 9. The comparison of the z-flow distributions for (x) = 0.0049 and (Q^) 
calculated with (solid) and without (dashed) the kinematic correction (|90|). 



32.6 GeV^ 



We should emphasize that we don't know the exact behavior of g^'^\b,x) for x > 10^^, 
where currently there is no data available. On the other hand, the term g^'^\b,x) in the 
parameterization (^) becomes negative for x > 5 ■ 10~^ which makes the numerical value of 
the resummed z-How (|73|) unphysical. Thus, the parameterization ( PHj ) must be used only 
at X < 10-2. 

The theoretical results in Fig. ^ were obtained using the kinematic correction to the 
asymptotic and resummed cross-sections at non-zero qt, which was discussed in Section 0. 
As can be seen from Fig. P, without this correction the agreement between the resummation 
calculation and the data is still good in the region ^t < 2 GeV < Q, where the resummation 
calculation is truly applicable. The kinematic correction improves the agreement between 
the resummed ^-flow and the data in the region gr > 2 GeV. In this region, the theoretical 
prediction without the kinematic correction significantly exceeds the data. As explained in 
Section |V|, this can be attributed to differences between the phase space of the perturbative 
and resummed pieces. In the case of the resummed cross-section the phase space may expand 
to much lower x and z than is allowed for the perturbative cross-section. Consequently, the 
large magnitude of PDFs and FFs at small x and z may spoil the cancellation between the 
resummed and asymptotic piece at large qx- This difference can be corrected by redefinition 



(pop of the variables x and z in the resummed and asymptotic pieces. 

One of the advantages of the resummed cross-sections is that, by their construction, 
they are less dependent on the choice of the renormalization and factorization scales of the 
problem, and on the end points Ci/b and C2Q in the integral of the perturbative Sudakov 
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FIG. 10. The dependence of the resumnied z-flow on the perturbative scale fi and the constants 
Ci,C2- The Figure corresponds to (x) = 0.0049, (Q^) = 32.6 GeV^. The curves correspond to: 
Ci = 2e~^, C2 = e"3/4,^ = Q (sohd); Ci = 2e-^, C2 = e'^/^/i = 2Q (dashed); Ci = 26"^, 
C2 = e-3/^^ = 0.5Q (dotted); Ci = 26"'^, C2 = 1,/x = Q (dot-dashed). 



factor (1^). The reason is that the scale variations in the perturbative part of ( [74D are 
compensated by the variation of the term g'^'^\b,Q) in the nonperturbative Sudakov factor 



In Fig. |10| we show the resummed 2;-flow for different choices of perturbative renormal- 
ization and factorization scale /i, varying between 0.5Q and 2Q, and for a different choice 
of the constants Ci and C2 (the "canonical choice" Ci = 2e~'^, C2 = 1). As expected, the 
resummed cross-section shows little variation with the changes oi fi, Ci, C2. 

We have also checked the stability of the resummed z-flow under variation of the momen- 
tum transfer scale Qo in the logarithm of g^'^\b,Q), and under variation of the parameter 
bmax separating the perturbative and nonperturbative dynamics in (ffil). Fig. |ll] shows the 
2-flows for variations of Qq and b^ax by factors of 1/2 and 2. As can be seen in Fig. |TI|, the 



variation of the resummed ^-fiow is small. The stability of the resummed ^-fiow under the 
variation of bmax, Qo,Ci, C2 indicates that the choice of these parameters has less influence 
on the shape of the resummed cross-section than the free parameters in g^^\b,x). It also 
illustrates the fact that the existing HERA data is relatively insensitive to the parameteriza- 
tion of the Q-dependent function g^'^\b, Q), which can be studied in a more detailed fashion 
once the data in a larger range of Q become available. 

Finally, in Fig. O we replot the results of our calculation presented in Fig. H as the cm. 



pseudorapidity distributions of the transverse energy flow (Et). This quantity is obtained 
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FIG. 11. Dependence of the resummed z-flow on Qq and bmax- 
to Qo = 1 GeV, bmax = 1 GeV"^ (thick sohd); Qq = I GeV, bmax 



The curves correspond 
= 0.5 GeV-i (dashed); 



Qo = l GeV, bmax = 2 GeV"^ (dotted); Qo = 0.5 GeV, bmax = 1 GeV"^ (thin sohd); Qo = 2 GeV, 
bmax = 1 GeV"^ (dot-dashed). The graph corresponds to (x) = 0.0049, (Q^) = 32.6 GeV^. 

by the transformation (0). The small-gr region, where the resummation formaUsm is valid, 
corresponds to large pseudorapidities. In this region, the agreement between our calculation 
and the data is good. At smaller pseudorapidities (larger qx), one sees the above-mentioned 
excess of the data over the perturbative NLO calculation. In the (Et) vs. r]cm plot, this 
excess is magnified because of the factor g|, in the transformation ([T2|). 
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FIG. 12. The 7*p cm. pseudorapidity distributions of the transverse energy flow in the current 



region. The data are from |11| 
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FIG. 13. The charged particle multiphcity at (x) = 0.0044, (Q^) = 35 GeV^ 
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B. Multiplicity of charged particle production 

The HI and ZEUS collaborations recently published the pseudorapidity distributions 
of charged particle multiplicity at small values of x and Q^ P]. The data presented was 
organized in the same bins of x and Q^ as those used in the HI analysis of the transverse 
energy flow |Tl|]. Also, the experiment E665 at Fermilab has provided extensive data on 



charged particle multiplicity covering the region of larger x [|10| . However, a comparison of 
these data with our calculation, analogous to the one described in the previous subsection, 
is complicated by several obstacles. 

In comparison to energy flows, charged particle multiplicity depends on the additional 
variable z which controls the fragmentation of a final-state parton into the observed hadrons. 
Therefore the pQCD hadropro duct ion cross-sections depend on parton fragmentation func- 
tions -D(^fe, /i), which at present are known well only in the region 0.1 < ^b < 0.8 p6 |. 



The non-perturbative Sudakov factor (^) for particle multiplicity can depend on z. Fur- 
thermore, in principle the non-perturbative Sudakov factor can be different for different types 
of initial and final states (see the discussion of Eq. (|6^) above). As a first approximation, in 
the analysis below we will ignore this difference. 

On the other hand, in the existing data on charged particle pseudorapidity or transverse 
momentum distributions |^|T^ the dependence on the final-state fragmentation variable is 
not separated from the other variables. Therefore more complicated fitting of the data is 
necessary to disentangle the qr and z dependences. 

More importantly, our calculation was made under the assumption that all the participat- 
ing particles, including the final-state hadrons, are massless. Because of this assumption, the 
production of soft final-state hadrons, with = is allowed. This contradicts the situation 
of the experiments at HERA and Fermilab, in which there is a non-zero minimal value of z 
determined by the finite mass of the observed hadron. It follows from the definition (^) of 
z and the formulas (^, |2^) for the initial and final hadron momenta in the 7*p cm. frame, 
that 

. = ^>^, (94) 

where 

p'^ = E"^ + p'Z and W^ = Q^(--l 



Using (p4D, one can show that the multiplicity distributions presented in [^ receive con- 
tributions from charged pions with z ~ 10~^ or lower. This means that for the experimental 
cuts used in the HI analysis, the multiplicity receives significant contribution from the region 
2; < 0.1 where the fragmentation functions are poorly known, and mass effects are important. 



In Fig. |T3| we present the resummed multiplicity in the region z > 0.1, where the uncer- 
tainties in our knowledge of the fragmentation functions and the mass effects are minimal. 
The results in Fig. 13 are given for x = 0.0044, Q'^ = 35 GeV^ and various bins of z. In the 



calculation, we used the CTEQ4M PDFs [^ and the fragmentation functions from |^ . We 
show the prediction using the z-independent non-perturbative Sudakov factor (|9ll) obtained 
from the study of the energy flows. We also assume that the charged particle production 
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rate is dominated by fragmentation of partons into charged pions, kaons and protons, and 
that the non-perturbative Sudakov factors for these types of particles are the same. In spite 
of many assumptions that go into this calculation, it will be interesting to compare it with 
the experimental data. 

VII. CONCLUSIONS 

In this paper, we have presented a formalism for all-order resummation of large logarithms 
arising in hadroproduction in the current region of deep-inelastic scattering, i.e. for large 
pseudorapidity of the final-state hadrons in the photon-proton cm. frame. We found that the 
formalism describes well the behavior of the transverse energy fiows measured at HERA |ll 



in the region of large hadronic cm. pseudorapidity r/ > 3.0. At smaller pseudorapidities, we 
found a deficit of the NLO rate compared to the existing data. Evidently, this is a signature 
of the importance of the NNLO corrections, which were not studied in this paper. 

The formalism presented here can also be directly applied to the study of hadron mul- 
tiplicities. In this case, however, additional care is required in treating the final-state frag- 
mentation of partons and the effects of the mass of the final-state hadron. In view of this, 
reanalysis of charged particle production multiplicity measured by HI, ZEUS and Fermilab- 
E665, with the goal to separate the z and g-r dependence, will be very useful to study the 
properties of the non-perturbative part of the Sudakov factor, and to improve the applica- 
bility of the perturbative calculation in the current region. 
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APPENDIX A. TRANSFORMATION FROM THE HADRON FRAME TO THE 

HERA LAB FRAME 

In this Appendix, we summarize the relationships between the relativistic invariants used 
in this paper, the hadron frame variables, and the particle momenta in the HERA lab frame. 

The definition of the HERA lab frame is that the proton {A) moves in the +z direction, 
with energy Ea-, and the incoming lepton moves in the —z direction with energy E. The 
momenta of the incident particles are 

P^ = (Ea,0,0,E^), (95) 

P = (E,0,0,-E) . (96) 
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The outgoing lepton has energy E' and scattering angle 9 relative to the —z direction. 
We define the x-axis of the HERA frame such that the outgoing lepton is in the a;2;-plane; 
that is, 



I'l" = (E', -E' sin 6, 0, -E' cos ( 



(97) 



The observed hadron {B) has energy Eb and scattering angle 9b with respect to the +z 
direction, and azimuthal angle 0b; thus its momentum is 



p^Qi = {Eb, Eb sin 9b cos 0^, Eb sin 9b sin 0^, Eb cos 9b) ■ 



{91 



The scalars x and Q^ are completely determined by measuring the energy and the scat- 
tering angle of the outgoing lepton, 



Q' 



X 



2EE'{l-cos9), 
EE'(l-cos( 



Ea[2E-E'{1 + cos9)] ■ 
The scalars z and q^ depend on the outgoing hadron and lepton as 
f3EB{l- cos 9b) 



(99) 
(100) 



Qt 



Q 

2Q^{l-{y-l)P^) 
P^{1- cos 9b) 



sm 



2 ^B 



+ sin 9b sin 9^ sin^ — 



where 



V 






(101) 
(102) 

(103) 



The angle and the boost parameter -0, which are used for the angular decomposition 
of the cross-sections in the hadron Breit frame, can be found as 



cos I 



Q 



2 VI -y?T 



cot\9B/2) ql 



/32 



Q' 



cosh-0 = 1. 

y 



(104) 
(105) 
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APPENDIX B. THE PERTURBATIVE CROSS-SECTION AND Z-FLOW 

DISTRIBUTION 



In this Appendix, we collect the formulas for the NLO parton level cross-sections 
dcTba/ {dxdzdQ'^dq^dcl)) . 

According to (|30|), the hadron level cross-section da ba/ {dxdzdQ'^ dq^d(j)) is related to the 
parton cross-sections da^a / {dxdzdQ"^ dq^dcj)) as 



dxdzdQ^dq^d(f) Tt -^^ ^^ "^^ ^^ " " dxdzdQ^dq^c 



At non-zero qr, the parton cross-section receives the contribution from the real emission 
diagrams (Fig. |e-f); it can be expressed as 



d&ba 



dxdzdQ'^dq^dcf) 
o-qFi a. 



AtiSsaQ'^ 71" 



4 fl ^Y]._, 



g2 \x 



J j a=l 



(106) 



with the same notations as in Section |T|. In this formula, 

J2f^'^\x,z,Q\q'^)A^ii;,<P) = 2S,,CpXz\[\(^ + iQ'-q'^r)+6QYi 



+ 2Q\2A2 + A4) + 2-^(Q2 + g2)^3 1 



(107) 



^/jf(x,^,Q^g2)^„(^,</)) 
0=1 

x{l — x)' 



2 , .2-^2 ..+2)+2Q\^-\-\] 
qj^ Kx'^z'^ xz J X z 



A, 



+ 4Q2(2A2 + A4) + 2-^(2(g2 + g|) - ?,)A^ 

qx xz 



(108) 



0=1 
2C^x(l - 5 



^/(f(x,^,Q^g^)A„(^,0) 



1/ Q^ 



g|. \x'^{l — zy 



+ (Q' - ^t) + 6Q' Ai 



In (Hi, 



qr 



2Q\2A^ + A4) + 2^(g2 + ^~2)^ 1 _ 



\-z 



(109) 



(110) 
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The index j corresponds to a quark (antiquark) of a type j, the index g corresponds to a 
gluon. 

From (pOD, it is possible to derive the perturbative 2;- flow distribution, 



dxdQ'^dq?pd(j) ^ Jz^i„ dxdzdQ'^ dqj^c 



r^%- E E e? /' ^F.IA{ia)z'x E A:^ {X, z, Q\ ql)AM, </>). (Ill) 

It depends on the same functions fba{x, z, Q^, q^), with the parton variable z determined by 
the (5-function in fO!Of 



^ ^ (112) 



{q^/Q^ ~l)x + l 
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